The Quantum Interference Effect Transistor 
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We propose a new type of molecular transistor, the Quantum Interference Effect Transistor 
(QuIET), based on tunable current suppression due to quantum interference. We show that any 
aromatic hydrocarbon ring has two-lead configurations for which current at small voltages is sup- 
pressed by destructive interference. A transistor can be created by providing phase relaxation or 
decoherence at a site on the ring. We propose several molecules which could tunably introduce 
the necessary dephasing or decoherence, as well as a proof of principle using a scanning tunneling 
microscope tip. Within the self-consistent Hartree-Fock approximation, the QuIET is shown to have 
characteristics strikingly similar to those of conventional field effect and bipolar junction transistors. 



I. INTRODUCTION 

Although there has been considerable scientific and 
commercial interest, a small molecular transistor has yet 
to be discovered. Needless to say, such a device is cru- 
cial to transferring existing technology to smaller length 
scales. Current industrial fabrication techniques have 
more or less exhausted the possibilities for purely clas- 
sical phenomena to solve this problem. We must there- 
fore turn to quantum mechanical effects in the search 
for smaller transistors. The solution, therefore, is also of 
fundamental interest to mesoscopic and molecular physi- 
cists: we find that the "small transistor" problem is an 
engaging way of posing the question, "What happens to 
multi-lead, many-body electronic systems when their size 
is such that quantum effects are important?" 

We propose a new type of device, the Quantum In- 
terference Effect Transistor (QuIET), capable of filling 
the role of the Field Effect Transistor (FET) and Bipolar 
Junction Transistor (BJT) at length scales <lnm. The 
QuIET consists simply of a hydrocarbon ring and a mech- 
anism to tunably introduce phase relaxation at a partic- 
ular site. Unlike the Single Electron Transistor (SET), 
the transistor behavior can occur over a large range of 
base voltages. Its / — V characteristic is a single, broad 
resonance, strikingly similar to those of macroscopic tran- 
sistors over a domain of several volts. This is to be taken 
in contrast to the SET's / — V, a series of many sharp 
peaks. 

The operating principle of the QuIET is that bias ap- 
plied to a third lead can modulate an otherwise complete 
conductance suppression across a hydrocarbon ring. This 
conductance suppression is a simple, single-particle effect 
of quantum mechanics. The modulation from a third lead 
can be achieved either through direct coupling to the 
molecule (Fig.^), as in the case of a scanning tunneling 
microscope (STM) tip, or via the introduction of an ap- 
propriate intermediary molecular complex, as shown in 
Fig.^. In either case, the resultant device is an excellent 
molecular-scale transistor. 

Section ^1 gives an outline of the extended Hubbard 
model Hamiltonian and Green function formalism we use 
to treat the problem. The multi-terminal current for- 
mula allows extraction of current. Section IlIII explains 




FIG. 1: Schematic diagrams of two types of QuIET. In each, 
base voltage modulates the coherent suppression of current 
between emitter (E) and collector (C) leads. In (a), base 
voltage controls the distance x between the benzene ring and 
base lead (B), for example an STM tip. This in turn controls 
the coupling of the ring to the base lead. In (b), a base 
complex is introduced between the ring and base lead. The 
electrostatic effect of the base lead's bias on this molecule 
alters its coupling to the benzene ring. 



the tunable coherent conductance suppression by which 
the QuIET works. Section IIVI focuses on the use of a 
acceptor-donor molecule as a scalable method of tuning 
interference. Section [V] presents numerical results which 
indicate the viability of this idea. Finally, Section IVII 
contains our conclusions. 



II. MODEL 

The Hamiltonian of the system can be written as the 
sum of three terms: H = Hm + Hi + Hmn- The first is 
the extended Hubbard model molecular Hamiltonian 

H = Y. edU^a^ E *v {dld.a + H.c.) +^ ^QiQ„ 
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where dia- annihilates an electron on atomic site i with 
spin a, Si are the site energies, and tij are the tunneling 
matrix elements. The final term of Eq. contains in- 
tersite and same-site Coulomb interactions, as well as the 
electrostatic effects of the leads. The interaction energies 



are modeled according to the Ohno parameterization 
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where Rij is the distance between sites i and j. Qi is an 
effective charge operator for atomic site i: 



1. 



(3) 



The second term represents the polarization charge on 
site i due to capacitive coupling with lead a. Here Cia 
is the capacitance between site i and lead a, chosen to 
correspond with the interaction energies of Eq. |(2Jl , and 
Vq is the voltage on lead a. e is the magnitude of the 
electron charge. 

The QuIET is intended for use at room temperature 
and above, a temperature range far beyond the regime 
in which lead-molecule or lead-lead correlations play an 
important role. As such, we have followed the method 
of Ref . i and treated electrostatic interactions between 
molecule and leads at the level of capacitance parame- 
ters. The electronic situation of the leads is thus com- 
pletely determined by the externally controlled voltages 
Va, along with the leads' temperatures and Fermi ener- 
gies. Each lead possesses a continuum of states, and their 
total Hamiltonian is 



H, 
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where ej. are the energies of the single-particle levels k 
in lead a, and Cka- are the annihilation operators for the 
states in the leads. 

Tunneling between molecule and leads is provided by 
the final term of the Hamiltonian: 



= EE(^'^-4-^'=- + H-c. 



(5) 
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Vik are the tunneling matrix elements for moving from a 
level k within lead a to the nearby site i. Coupling of the 
leads to the ring via inert molecular chains, as may be 
desirable for fabrication purposes, can be included in the 
effective Vik , as can the effect of the substituents used to 
bond the leads to the molecule. 



A system whose Hamiltonian includes such terms as 
Eqs. 10} and |(SJ) requires an infinite-dimensional Pock 
space, but we are concerned mainly with the behavior 
of the discrete molecule suspended between the leads. 
We therefore adopt a Green function approach, in which 
Dyson's Equation gives the Green function of the full 
system 



GiE)=[G-\E)-^{E)]- 



(6) 



where Gm is the Green function of the isolated molecular 
system. With the use of an appropriate self-energy E, 
Equation ((BJ is true both for the retarded Green function 
G"' as well as for its 2 x 2 Keldysh counterpart. 
The retarded self-energy due to the leads is 



(7) 



where the energy widths are given by Fermi's Golden 
Rule 



r,(i?) = 27r^^|K,.p<5(£;-4), 



(8) 



We take the broad-band limit of Eq. ||SJ) and treat 
each of the F^ as a constant parameter characterizing 
the lead-site coupling. The only effect of Hi and Htun 
in this limit is to shift the poles of the Green function 
into the complex plane. This causes the density of states 
p{E) = — ilmTrG"'(£') to change from a discrete spec- 
trum of delta functions to a continuous, width-broadened 
function. Due to the open nature of the system, electrons 
can occupy all energies. 

The retarded Green function gained via Eq. © con- 
tains all information regarding the dynamics of the sys- 
tem. In particular, the current in lead f3 is given by the 
familiar multi-terminal current formula 0: 



dET^0{E)[UE)~fp{E)], (9) 



where fa is the Fermi function for lead a. The transmis- 
sion probability is 



Taf,{E) = rart\G:,{E)\ 



(10) 



Here a{b) is the site with hopping to lead a(/3). We note 
that Eq. ^ is an exact result of the Keldysh formalism 
in cases, like ours, consisting only of elastic processes. 

In order to arrive at G^{E), we must consider electron- 
electron interactions. Here we do so via the well known 
self-consistent Hartree-Fock method. The Hamiltonian 
is replaced by its mean-field approximation 
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In this approximation, the retarded Green function is 
G^^{E) = {E-H^^^^O+)-\ 

Equation gives the mean-field Hamiltonian as a 
function of the diagonal and off-diagonal equal-time cor- 
relation functions {d\^djcr). To complete a self-consistent 
loop, we require an expression for these quantities in 
terms of Green functions. They are given in the Keldysh 
formalism by the equal-time limit of the "<" Green func- 
tion 

/OO 7 
^G<^^.,(c.)e— (*- 

From Dyson's Equation ijHJl it follows 01 that 

G'< = G"'I]<G"'t + (1 + G''I]'^)G< (1 + S'-tC^). (13) 

The second term is a purely equilibrium property, and 
can be related to the total charge on the molecule when 
the three lead biases are equal. The "<" self-energy is 
given by 

^<,{u)^lTJo.{L0)5ab. (14) 

The desired relation between the equal-time correlation 
functions and C, 

/OO 7 

(15) 

is now readily computed, and the self-consistent loop is 
complete. 

III. TUNABLE CONDUCTANCE 
SUPPRESSION 

In the two- lead device shown in Fig.[21i, a single carrier 
(electron or hole) is injected into the ring by the emitter 
and exits via the collector some time later. In the path 
integral formulation of quantum mechanics, it traverses 
all paths around the ring allowed by the connectivity of 
the system during this process. 

We operate the QuIET in the regime where there is 
little charge transfer between it and the leads. In the 
linear response, the carrier has momentum equal to the 
Fermi momentum of the ring kp — j;^, where a — 1.397A 
is the intersite spacing of benzene. Clearly, then, the 
phase difference between paths ® and ® is tt, and they 
cancel exactly. Similarly, all of the paths through the 
ring from emitter to collector exactly cancel in a pairwise 




FIG. 2: (a) Two-lead experiment to measure the conduc- 
tance of benzene when the leads are in the meta configuration. 
Shown are the two most direct paths a carrier can take from 
the emitter lead to the collector. These two cancel exactly, 
as do all other paths with the same endpoints in a similar 
pairwise fashion, (b) Example of a new path allowed when 
a base complex or lead is included. Such paths are not can- 
celed, and so contribute to the total current between emitter 
and collector. 



fashion. Therefore, transport of carriers is forbidden in 
linear response. 

It is a consequence of Luttinger's Theorem p] that 
this coherent suppression of current persists into the in- 
teracting regime, as demonstrated in Fig.|2K. The trans- 
mission is calculated by the self-consistent Hartree-Fock 
model outlined in Section ^1 In this figure the base cou- 
pling Fb = 0, and so transmission at the Fermi energy 
is wholly suppressed by coherence. Figure l^) shows how 
this result changes as F^ is increased. Current is al- 
lowed to flow due to two new phenomena. The first is 
that new paths, such as the one shown in Fig. are 
added. These paths have no particular phase relation- 
ship to other paths, and so support transmission. Addi- 
tionally, nonzero Tb denotes decoherence, and therefore 
a departure from the perfect coherent current suppres- 
sion of Fs = 0. This is the basic operating principle of 
the QuIET: coherent current suppression can be tunably 
broken by the introduction of decoherence and dephas- 
ing from a base complex or third lead. In fact. Figure 
13: shows that the transmission varies nearly linearly for 
Fs < t. 

Clearly, the transistor behavior based on the mech- 
anism outlined above is requisite upon an assumption 
that the device be operated well within the gap of ben- 
zene. The numerical simulations discussed in Section 
IVl indicate that the best transistor results are found for 
collector-emitter bias ^ 1 — 2V. Another, related, consid- 
eration is that in equilibrium, charge transfer between the 
molecule and leads not play an important role. For this 
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FIG. 3: Total transmission probability of the device shown in Fig. Tec = Tec + Teb+Tbc ' ^^^'^^ second term is due 
to a restriction that the base lead be an infinite-impedance ideal voltage probe 0- Here, Te = and Vc = ■2t, but the 
choice of energy widths does not affect the qualitative results. In (a), Tb ~ 0. In (b), the lowest curve shows Fs ~ and each 
successive curve shows increasing F_b in increments of .It. (c) shows the transmission at the Fermi energy as a function of Fs. 




FIG. 4: Emitter-collector lead configurations possible in a 
QuIET based on [18]-annulene. The bold lines represent the 
positioning of the two leads. Each of the four arrangements 
has a different phase difference associated with it: (a) n, (b) 
Stt, (c) 5-k, and (d) 7n. 



to be true, the work function of the metalhc leads must 
be comparable to the chemical potential of benzene. For- 
tunately, this is the case with many bulk metals, among 
them palladium, iridium, platinum, and gold 

While in this work we have chosen to focus on benzene, 
the QuIET mechanism applies to any aromatic annulene 
with leads positioned so the two most direct paths have 
a phase difference of tt. Furthermore, larger molecules 
have other possible lead configurations, based on phase 
differences of Stt, Stt, etc. Figure^jshows the lead config- 
urations for a QuIET based on [18]-annulene. Of course, 
benzene provides the smallest of all possible QuIETs. 

The position of the third lead affects the degree to 
which destructive interference is suppressed. For ben- 
zene, the most effective location for a third lead is shown 
in Fig.[S^. The base may also be placed at the site imme- 
diately between the emitter and collector leads, as shown 
in Fig. The QuIET operates in this configuration as 
well, although since base coupling to the current carri- 
ers is less, the transistor effect is somewhat suppressed. 
The third, three-fold symmetric configuration of leads 
(Fig. O;) completely decouples base from current carri- 
ers within benzene. Because of this, the base cannot 
be used to provide the decoherence or dephasing neces- 



(a) (b) (c) 




FIG. 5: The three different arrangements for a third lead on 
benzene when the emitter and collector leads (bold) are in 
the meta configuration, (a) is the choice which couples most 
strongly to the conducting orbitals of benzene, (b) is a second 
case which allows for QuIET operation. The third possibility, 
(c), decouples entirely from the conducting molecular orbitals 
by symmetry. A third lead in this configuration cannot break 
the coherent current suppression at all, and so the molecule 
does not function as a QuIET. 

sary to QuIET operation in this configuration. For each 
aromatic hydrocarbon, exactly one three-fold symmetric 
lead configuration exists and yields no transistor behav- 
ior. 



IV. BASE COMPLEX 

With the tunable current suppression outlined in the 
previous section, the working principle of a transistor has 
become apparent. The most straightforward method of 
varying F^ is to change the distance of the third lead 
from the hydrocarbon molecule, as via an STM tip. It is 
also possible, however, to interpose an additional molec- 
ular complex between annulene and the base lead. If the 
effect of biasing the base lead is to increase the trans- 
parency of the molecular complex, dephasing and deco- 
herence result. Such a QuIET is depicted schematically 
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FIG. 6: Specific example of a QuIET: a ditiiiol derivative 
of phenylene-TCNQ. E and C are the emitter and collector 
leads. The base lead can be either weakly coupled via TCNQ, 
as for a BJT analogue, or capacitively coupled if FET-like 
behavior is desired. 



in Fig.^. 

A simple way to engineer the base complex is to in- 
clude two spatially separated orbitals, with a detuning 
from each other A and hopping to each other t^. If 
A ^ is for zero bias on the base lead, paths through 
the base complex are highly suppressed. Since the or- 
bitals are spatially separate, however, one electrostati- 
cally couples more strongly to the base lead. Thus, the 
voltage on that lead can be used to control A, and hence 
change the transparency of the base complex. The re- 
sult is the broad conductance peak characteristic of two 
discrete levels coming into resonance. 

One way to achieve this structure is simply to link 
the annulene via an inert bridge to a donor or accep- 
tor molecule, e. g. the dithiol derivative of phenylene- 
TNCQ shown in Fig. O In this case, the effect of a base 
bias of appropriate sign is to bring the lowest unoccupied 
molecular orbital (LUMO) of TCNQ into resonance with 
the neighboring tt orbital of phenylene. The extra paths 
allowed within and through the TCNQ break the coher- 
ent suppression of current within the benzene ring. Any 
donor or acceptor molecule with which phenylene bonds 
in this way will yield similar results, although a donor 
used in this way will allow emitter-collector current to 
flow for negative, rather than positive, base voltages. 

The class of donor-acceptor rectifiers proposed by Avi- 
ram and Ratner |^ have this two-level structure as well. 
These molecules possess an orbital of increased electron 
affinity (the acceptor) and one of increased ionization po- 
tential (the donor), as well as perhaps an inert bridge 
orbital. By the influence of a base lead, the donor high- 
est occupied molecular orbital (HOMO) and acceptor 
LUMO can be brought into resonance, and the rectifier 
be made to serve as an effective QuIET base complex. 
In addition to this effect, several many-body mechanisms 
have been proposed whereby such molecules could exhibit 
asymmetric current flow For our purposes, however, 
it is sufficient to note that rectiflcation in general is in- 
dicative of a tunable transparency to charge carriers. 



TABLE L Approximate ranges of values available for molec- 
ular diodes which could be used as QuIET base complexes. 



Parameter 


Range 


Donor HOMO energy 


[-8,-5]eV 


Acceptor LUMO energy 


[-3,-l]eV 


tB 


[.01,l]eV 


t' 


[.01,l]eV 


Distance between donor and acceptor 


[6,10]A 



The variety within this well studied family of molecules 
(e. g., Refs. 0, H [H E d ) lends a great deal of 
versatility to the QuIET. Different choices of donor, ac- 
ceptor, and bridge molecules allow QuIETs to be fabri- 
cated for specific applications. Table^gives approximate 
ranges of parameters available for the use of molecular 
diodes as base complexes. 



V. NUMERICAL RESULTS AND DISCUSSION 

We turn now to discussion of numerical results, based 
on the self-consistent Hartree-Fock model presented in 
Sectionin For the purposes of illustration, the numerical 
results are based on the dithiol derivative of phenylene- 
TCNQ shown in Fig. Inland three bulk gold leads. Similar 
results can be obtained for many such QuIETs. 

For the purposes of electrostatic interactions, each lead 
is considered to be immediately adjacent to its neighbor- 
ing site. The electrostatic interactions between TCNQ 
and phenylene are modeled as those between two con- 
ducting spheres of appropriate size, while within pheny- 
lene the Ohno parameterization ((SJ is used. We neglect 
reduction of the Uij due to screening from a third site, an 
effect which we find to be of the order of 1%. In keeping 
with our assumption that the leads do not significantly 
perturb the molecular system, modifications to the inter- 
action parameters due to the presence of the leads are 
also neglected. 

A typical I — V diagram for this QuIET is shown in 
Fig. [7| demonstrating that the QuIET is quite reminis- 
cent in operation to classical transistors. The currents 
in the emitter and collector leads exhibit a broad reso- 
nance as the base voltage is increased. Furthermore, for 
nonzero Tb, the device amplifies the current in the base 
lead, providing emulation of the classical BJT. Since a 
donor, TCNQ, was used in the base complex, the device 
performs in a manner analogous to an NPN transistor, 
with base current flowing into the molecule. Use of an 
acceptor instead, for example TTF, yields PNP-like be- 
havior. For capacitive coupling of the third lead, on the 
other hand, a FET-like I — V characteristic is obtained. 

We interpret this transistor behavior as due to the co- 
herence mechanism discussed in Section IIIII If hopping 
between the benzene ring and the base complex is set to 
zero, we find that full coherent current suppression is re- 
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FIG. 7: Typical I — V characteristic of two QuIETs, show- 
ing emitter current out of the molecule vs. voltage on the base 
lead for three different biases applied to the collector, with the 
emitter ground. The calculation is done for the molecule in 
Fig.|Slat room temperature. T = ^^+^0 Sives the sequential 
tunneling rate through the device. Here, Te = Tc = 2.4eV. 
In (a), Fb = .0024eV, and the QuIET amplifies current, sim- 
ilar to a BJT. The dotted curve is the base current into the 
molecule for the case of Vce ~ .24V, multiplied by 10 for 
clarity. Base currents for other collector voltages are similar, 
(b) shows the FET limit of the QuIET, with Fs = 0. 
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FIG. 8: Ie vs. Vbc for different values of Fb = Fc = F. Each 
curve has been scaled so that the leakage current at zero bias 
is 1. Fb = 0, Vce = .24V, and positive current is out of 
the molecule. Reducing F suppresses the leakage current, but 
does not strongly affect the QuIET's behavior on-resonance, 
thus greatly enhancing its overall function. 



stored and almost no current flows between the emitter 
and collector leads. Furthermore, the current step ef- 
fect persists for arbitrarily small Fb, which is consistent 
with our interpretation that transport through the non- 
canceling paths is enhanced by the electrostatic effect of 
the third lead. 

Estimation of the F^'s remains an open question in the 
field of molecular electronics. Similar quantities are of- 
ten estimated to be ^ .5eV 0, by the method of 
Ref. [la. w hereas values as high as leV have been sug- 
gested [131 • For the emitter and collector broadenings, 
we have taken values somewhat higher than the norm for 
the sake of numerical convergence. Fortunately, nothing 
in the arguments of Sections IIIII or IIVI depends strongly 
on the choice of these quantities. So long as they are pos- 
itive, their magnitude merely determines the scale of the 
current, and has negligible effect on the overall operating 
principle and transistor behavior of the QuIET. 

Figure |S1 demonstrates the effect of varying the cou- 
pling of the benzene ring to emitter and collector leads. 
As is to be expected, the internal structure of the res- 
onance, due to electrons moving to screen the chang- 
ing field from the base, sharpens. Furthermore, at low 
base voltages, less than about .75V, the base complex is 
strongly off resonance, and the scale of the leakage cur- 
rent through the annulene is set mainly by F = -f^^f^, 
an expected result for systems like the QuIET, which are 
not near a charge fluctuation resonance . 

For higher base voltages, the base complex begins to 
play an important role. Thus, as described by Ref. |l5j |. 
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r alone no longer determines the scale of the current. 
Instead, we find that the rate-limiting process is travel 
through the base complex, and varying Te and Tc has 
little effect. Thus, smaller emitter and collector energy 
widths enhance the QuIET's current contrast dramati- 
cally. 

VI. CONCLUSIONS 

We have presented a novel idea for a small single- 
molecule transistor. The device is based on the coherent 
suppression of current through an aromatic hydrocarbon 
ring, and the control of that effect by increasing the con- 
tribution of paths outside the ring. Two simple methods 
of creating that control, either through an STM tip or 



through a small base complex, were discussed, and nu- 
merical results presented. 

In contrast to SET-style devices, the QuIET is pre- 
dicted to be an extremely versatile, scalable device. Since 
it is chemically constructed, it should be possible to fabri- 
cate as many identical QuIETs as desired. Furthermore, 
QuIETs can be designed which mimic the functionality 
of all major classes of macroscale transistors: FET, NPN, 
and PNP. 
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